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Resolving Cosmic Gamma Ray Anomalies with Dark Matter Decaying Now 
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Dark matter particles need not be completely stable, and in fact they may be decaying now. We 
consider this possibility in the frameworks of universal extra dimensions and supersymmetry with 
very late decays of WIMPs to Kaluza-Klein gravitons and gravitinos. The diffuse photon background 
is a sensitive probe, even for lifetimes far greater than the age of the Universe. Remarkably, both 
the energy spectrum and flux of the observed MeV 7-ray excess may be simultaneously explained by 
decaying dark matter with MeV mass splittings. Future observations of continuum and line photon 
fluxes will test this explanation and may provide novel constraints on cosmological parameters. 
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PACS numbers: 95.35.+d, ll.10.Kk, 12.60.-i, 98.80.Cq 

The abundance of dark matter is now well known from 
observations of supernovae, galaxies and galactic clusters, 
and the cosmic mocriwave background (CMB) but 
its identity remains elusive. Weakly-interacting massive 
particles (WIMPs) with weak-scale masses ~ 0.1 — 1 TeV 
are attractive dark matter candidates. The number of 
WIMPs in the Universe is fixed at freeze-out when they 
decouple from the known particles about 1 ns after the 
Big Bang. Assuming they are absolutely stable, these 
WIMPs survive to the present day, and their number 
density is naturally in the right range to be dark matter. 
The standard signatures of WIMPs include, for example, 
elastic scattering off nucleons in underground laborato- 
ries, products from WIMP annihilation in the galaxy, and 
missing energy signals at colliders 0. 

The stability of WIMPs is, however, not required to 
preserve the key virtues of the WIMP scenario. In fact, 
in supersymmetry (SUSY) and other widely-studied sce- 
narios, it is just as natural for WIMPs to decay after 
freeze-out to other stable particles with similar masses, 
which automatically inherit the right relic density to be 
dark matter 0]. If the resulting dark matter interacts 
only gravitationally, the WIMP decay is very late, in 
some cases leading to interesting effects in structure for- 
mation [4| and other cosmological observables. Of course, 
the WIMP lifetime depends on Am, the mass splitting 
between the WIMP and its decay product. For high de- 
generacies, the WIMP lifetime may be of the order of or 
greater than the age of the Universe to ~ 4.3 x 10 17 s, 
leading to the tantalizing possibility that dark matter is 
decaying now. 

For very long WIMP lifetimes, the diffuse photon back- 
ground is a promising probe 0, [B| . Particularly interest- 
ing is the (extragalactic) cosmic gamma ray background 
(CGB) shown in Fig. Q] Although smooth, the CGB 
must be explained by multiple sources. For E 1 < 1 MeV 
and Ej > 10 MeV, the CGB is reasonably well-modeled 
by thermal emission from obscured active gala ctic nu- 
clei (AGN) [9( and beamed AGN, or blazars [lfj, respec- 
tively. However, in the range 1 MeV < E 1 < 5 MeV, no 
astrophysical source can account for the observed CGB. 
Blazars are observed to have a spectral cut-off ~ 10 MeV, 
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FIG. 1: The CGB measured by HEAO-1 Q (circles), COMP- 
TEL (squares), and EGRET (triangles), along with the 
known astrophysical sources: AGN (long-dash), SNIa (dot- 
ted), and blazars (short-dash, and dot-dashed extrapolation). 



and also only a few objects have been detected below this 
v\ ; a maximal upper limit fl3| on the blazar 
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contribution for E~ f < 10 MeV is shown in Fig. [TJ Dif- 
fuse 7-rays from Type la supernovae (SNIa) contribute 
below ~ 5 MeV, but the most recent astronomical data 
show that they also cannot account for the entire spec- 
trum [3, EH ; previous calculations suggested that SNIa 
are the dominant source of 7-rays at MeV energies [la ]. 

In this paper, we study the contribution to the CGB 
from dark matter decaying now. We consider simple 
models with extra dimensions or SUSY in which WIMP 
decays are highly suppressed by both the weakness of 
gravity and small mass splittings and are dependent on 
a single parameter, Am. We find that the CGB is an 
extremely sensitive probe, even for lifetimes r 3> to- I n_ 
triguingly, we also find that both the energy spectrum 
and the flux of the gamma ray excess described above are 
naturally explained in these scenarios with Am ~ MeV. 

As our primary example we consider minimal univer- 
sal extra dimensions (mUED) fl7| . one of the simplest 
imaginable models with extra dimensions. In mUED all 



2 



particles propagate in one extra dimension compactified 
on a circle, and the theory is completely specified by m/,, 
the Higgs boson mass, and i?, the compactification ra- 
dius. (In detail, there is also a weak, logarithmic depen- 
dence on the cutoff scale A [l8|. We present results for 
AR = 20.) Every particle has a Kaluza-Klein (KK) part- 
ner at every mass level ~ m/R, m — 1,2,..., and the 
lightest KK particle (LKP) is a dark matter candidate, 
with its stability guaranteed by a discrete parity. 

Astrophysical and particle physics constraints limit 
mUED parameters to regions of (i? _1 ,m^) parameter 
space where the two lightest KK particles are the KK 
hypercharge gauge boson B 1 , and the KK graviton G 1 , 
with mass splitting Am < O(GeV) 19]. This extreme 
degeneracy, along with the fact that KK gravitons inter- 
act only gravitationally, leads to long NLKP lifetimes 
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where Mp ~ 2.4 x 10 18 GeV is the reduced Planck scale, 
6w is the weak mixing angle, b = 10/3 for B 1 — ► G 1 7, 
and 6=2 for G 1 — + B 1 ^ [20|. Note that r depends only 
on the single parameter Am. For 795 GeV < R^ 1 < 
809 GeV and 180 GeV < m h < 215 GeV, the model is 
not only viable, but the B 1 thermal relic abundance is 
consistent with that required for dark matter [2l[ and 
Am < 30 MeV, leading to lifetimes r(S x -> G 1 ^) > t . 

We will also consider supersymmetric models, where 
small mass splittings are also possible, since the gravitino 
mass is a completely free parameter. If the two light- 
est supersymmetric particles are a Bino-like neutralino B 
and the gravitino G, the heavier particle's decay width 
is again given by Eq. (Q]), but with b = 2 for B — ► G7, 
and b = 1 for G — > B7. As in mUED, r depends only on 
Am, and Am < 30 MeV yields lifetimes greater than to- 

The present photon flux from two-body decays is 
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where N(t) = N in e~ t/T and N in is the number of WIMPs 
at freeze-out, Vq is the present volume of the Universe, 
a is the cosmological scale factor with a(to) = 1, and 
e~ ( = Am is the energy of the produced photons. Pho- 
tons from two-body decays are observable in the diffuse 
photon background only if the decay takes place in the 
late Universe, when matter or vacuum energy dominates. 
In this case, Eq. $Z\) may be written as 
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where P(a) = t is the solution to (da/dt)/a = H(a) = 
HaVn M a- 3 + Que a- 3 ( 1 +™) with P(0) = 0, and Q M 
and f^DE are the matter and dark energy densities. If 



dark energy is a cosmological constant A with w = — 1 
_ 2 In [(v^aF +V^m + Oa« 3 ) A/H^] 
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The flux has a maximum at £l» 
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where [/(x) = (x + 1 - V3x + \)/{x - 1). 

The energy spectrum is easy to understand for very 
long and very short decay times. For r <C to, 
HqT 2 Hbe <C 1, and the flux grows due to the deceler- 
ated expansion of the Universe as d<fr/dE 7 oc E 1 / 2 until 
it reaches its maximum at E™ ax ~ £ 7 (fi M FT^T 2 /4) 1 / 3 . 
Above this energy, the flux is suppressed exponentially 
by the decreasing number of decaying particles 0] • 

On the other hand, if r 3> t a , HqT 2 ^!^ ^ U an d 
the flux grows as d^/dE 1 cx E 1 ^ 2 only for photons that 
originated in the matter-dominated epoch. For decays 
in the vacuum-dominated Universe, the flux decreases 
asymptotically as d^/dE^ oc _E( 1 + 3tu )/ 2 due to the accel- 
erated expansion. The flux reaches its maximal value at 
£7 ~ e 7 [-0 m /((1 + 3w)ftDE)r 1/(3M,) where photons 
were produced at matter- vacuum equality. Note that this 
value and the spectrum shape depend on the properties 
of the dark energy. Assuming Qm — 0.25, JIde = 0.75, 
w = — 1, and h = 0.7, and that these particles make up 
all of non-baryonic dark matter, so that 



1.0 x 10 _y cm" 



TeV 



in 



n 



NBDM 



0.2 



, (5) 



we find that the maximal flux is 

— (ff") = 1.33 x 10~ 3 cm" 2 s- 1 sr" 1 MeV" 1 
dE-y 7 
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Fig. [5] shows example contributions to the CGB from 
decaying dark matter in mUED and SUSY. The mass 
splittings have been chosen to produce maximal fluxes 
at E~ t ~ MeV. These frameworks are, however, highly 
constrained: once Am is chosen, r and the flux are es- 
sentially fixed. It is thus remarkable that the predicted 
flux is in the observable, but not excluded, range and 
may explain the current excess above known sources. 

To explore this intriguing fact further, we relax model- 
dependent constraints and consider r and Am to be in- 
dependent parameters in Fig. [3J The labeled curves give 
the points in (r, Am) parameter space where, for the 
WIMP masses indicated and assuming Eq. ([5]), the max- 
imal flux from decaying dark matter matches the flux of 
the observed photon background in the keV to 100 GeV 
range 0]. For a given WIMP mass, all points above the 
corresponding curve predict peak fluxes above the ob- 
served diffuse photon background and so are excluded. 

The shaded band in Fig.[3]is the region where the max- 
imal flux falls in the unaccounted for range of 1-5 MeV. 
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FIG. 2: Data for the CGB in the range of the MeV excess, 
along with predicted contributions from extragalactic dark 
matter decay. The curves are for B 1 — > G lr y in mUED with 
lifetime r = 10 3 to and m B i = 800 GeV (solid) and B — > G7 
in SUSY with lifetime r = 5 x 10 3 to and = 80 GeV 
(dashed). We have assumed Qnbdm = 0.2 and smeared all 
spectra with energy resolution AE/E = 10%, characteristic 
of experiments such as COMPTEL. The dot-dashed curve is 
the upper limit to the blazar spectrum, as in Fig. Q. 
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does not track Am, as the peak energy is significantly 
redshifted. For example, for a WIMP with mass 80 GeV, 
t - 10 12 s and Am ~ MeV, £™ ax ~ keV. The over- 
lap of this band with the labeled contours is where the 
observed excess may be explained through WIMP de- 
cays. We see that it requires 10 20 s < r < 10 22 s and 
1 MeV < Am < 10 MeV. These two properties may 
be simultaneously realized by two-body gravitational de- 
cays: the diagonal line shows the relation between r and 
Am given in Eq. |T]) for B 1 — > G 1 7, and we see that this 
line passes through the overlap region! Similar conclu- 
sions apply for all other decay models discussed above. 

These considerations of the diffuse photon background 
also have implications for the underlying models. For 
mUED, Am = 2.7 - 3.2 MeV and r = 4 - 7 x 10 20 s can 
explain the MeV excess in the CGB. This preferred region 
is realized for the decay B 1 -> G 1 7 for R- 1 w 808 GeV. 
(See Fig. 0J) Lower predicts larger Am and shorter 
lifetimes and is excluded. The MeV excess may also be 
realized for G 1 — * B 1 -/ for w 810.5 GeV, though in 
this case the G 1 must be produced non-thermally to have 
the required dark matter abundance [2^, • 

So far we have concentrated on the cosmic, or extra- 
galactic, photon flux, which is dependent only on cosmo- 
logical parameters. The Galactic photon flux depends 
on halo parameters and so is less robust, but it has 
the potential to be a striking signature, since these pho- 
tons are not redshifted and so will appear as lines with 
E~ ( = Am. INTEGRAL has searched for photon lines 
within 13° from the Galactic center 23j. For lines with 
energy E ~ MeV and width AE ~ 10 keV, INTEGRAL'S 



FIG. 3: Model-independent analysis of decaying dark mat- 
ter in the (r, Am) plane. In the shaded region, the result- 
ing extragalactic photon flux peaks in the MeV excess range 
1 MeV < £™ ax < 5 MeV. On the contours labeled with 
WIMP masses, the maximal extragalactic flux matches the 
extragalactic flux observed by COMPTEL; points above these 
contours are excluded. The diagonal line is the predicted re- 
lation between r and Am in mUED. On the dashed line, the 
predicted Galactic flux matches INTEGRAL'S sensitivity of 
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energy resolution at these energies, INTEGRAL'S sensi- 
tivity is <& ~ 10 -4 cm -2 s _1 . The Galactic flux from 
decaying dark matter saturates this limit along the ver- 
tical line in Fig. 02 assuming m x = 800 GeV. This flux is 
subject to halo uncertainties; we have assumed the halo 
density profiles of Ref . 24| , which give a conservative up- 
per limit on the flux within the field of view. Remarkably, 
however, we see that the vertical line also passes through 
the overlap region discussed above. If the MeV CGB 
anomaly is explained by decaying dark matter, then, the 
Galactic flux is also observable, and future searches for 
photon lines will stringently test this scenario. 

In conclusion, well-motivated frameworks support the 
possibility that dark matter may be decaying now. We 
have shown that the diffuse photon spectrum is a sensi- 
tive probe of this possibility, even for lifetimes r ^> to. 
This is the leading probe of these scenarios. Current 
bounds from the CMB [25| and reionization [2(| do not 
exclude this scenario, but they may also provide comple- 
mentary probes in the future. We have also shown that 
dark matter with mass splittings Am ~ MeV and life- 
times r ~ 10 3 — 10 4 Gyr can explain the current excess of 
observations above astrophysical sources at E 1 ~ MeV. 
Such lifetimes are unusually long, but it is remarkable 
that these lifetimes and mass splittings are simultane- 
ously realized in simple models with extra dimensional 
or supersymmetric WIMPs decaying to KK gravitons 
and gravitinos. Future experiments, such as ACT [27l |. 
with large apertures and expected energy resolutions of 
AE/E = 1%, may exclude or confirm this explanation of 
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FIG. 4: Phase diagram of mUED. The top and bottom shaded 
regions are excluded for the reasons indicated [ljj. In the 
yellow (light) shaded region, the B 1 thermal relic density is 
in the 2a preferred region for non-baryonic dark matter [2ll | . 
In the vertical band on the left (right) the decay B 1 — > G lr y 
(G 1 — > B 1 7) can explain the MeV diffuse photon excess. 

the MeV excess through both continuum and line signals. 
Finally, we note that if dark matter is in fact decaying 
now, the diffuse photon signal is also sensitive to the 
recent expansion history of the Universe. For example, 
as we have seen, the location of the spectrum peak is a 
function of £Im / ^de and w. The CGB may therefore, in 
principle, provide novel constraints on dark energy prop- 
erties and other cosmological parameters. 
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